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Abstract The Interplanetary Magnetic Field shows complex spatial and tem- 
poral variations. Single spacecraft measurements reveal only a one dimensional 
section of this rich four dimensional phenomenon. Multi-point measurements 
of the four Cluster spacecraft provide a unique tool to study the spatiotempo- 
ral structure of the field. Using Cluster data we determined three dimensional 
correlation functions of the fluctuations. By means of the correlation function 
one can describe and measure field variations. Our results can be used to verify 
theoretical predictions, to understand the formation and nature of solar wind 
turbulence. We found that the correlation length varies over almost six orders of 
magnitude. The IMF turbulence shows significant anisotropy with two distinct 
populations. In certain time intervals the ratio of the three axes of the correlation 
ellipse is 1/2.2/6 while in the remaining time we found extremely high correlation 
along one axis. We found favoured directions in the orientation of the correlation 
ellipsoids. 

Keywords: Magnetic fields, Interplanetary; Solar Wind, Disturbances; Turbu- 
lence 



1. Introduction and motivation 

The Interplanetary Magnetic Field (IMF) shows complex spatial and temporal 
variations which have been of much interest for a long time now. The phe- 
nomenon is often referred to as the turbulence of the solar wind. It is an 
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important factor in the energy and momentum transport of the plasma, in- 
fluences the formation and dynamics of shock waves and other discontinuities. 
The turbulence of the solar wind has its influence even on the propagation of cos- 
mic ray particles ( [Bieber, Wanner, and Matthaeus, 1996[ ) or on the distribution 
and mean free path of pickup ions dNemeth, Erdos, and Balogh, 2000] ). Every 
detailed model of the solar wind must account for these strong fluctuations. 

The magnetic variations of the solar wind are supposed to be the super- 
position of two types of fluctuations: the so called "2D-turbulence" and slab 
variations. In the model of slab disturbances the fluctuation wave vector is 
aligned with the mean magnetic field. This model has a physical motivation 
in terms of Alfven waves propagating along the mean field. In the 2D model 
however the turbulence wave vectors are predominantly in the plane perpen- 



dicular to the mean magnetic field (Matthaeus, Goldstein, and Roberts, 1990 



Bieber, Wanner, and Matthaeus, 1996 1. An interpretation of the 2D model using 



resonant three wave interactions was provided by Shebalin, Matthaeus, and Montgomery] 
(jl983p . More recently an alternative explanation was provided based on plasma 
mixing ( Nemeth et ai, 2006 1. 

The solar wind turbulence is known to be anisotropic and the contour plot 
of the two-dimensional correlation of the fluctuations as a function of distance 
parallel and perpendicular to the mean magnetic field shows a "Maltese cross" 



shape (Matthaeus, Goldstein, and Roberts, 1990 Matthaeus et ai, 2005). The 



Ulysses mission made it possible to study the turbulence in the high latitude 
solar wind as well ( [Bigazzi et ai, 2006 ). The anisotropy is present in the fast 
and slow solar wind alike ( Dasso et ai, 2005 ) . 

Single spacecraft measurements ( [Matthaeus, Goldstein, and Roberts, 1990 Bieber, Wanner, and Matthaei 
Dasso et ai, 2005 Bigazzi et ai, 2006) however only reveal a one dimensional 



section of this rich 3-1-1 dimensional phenomenon, which means that the spatial 
structure can only be estimated using strong assumptions. To avoid this limita- 
tion spacecraft couples ( Richardson and Paularena, 2001 ) were used at first. The 
multi-point measurements of the four Cluster spacecraft provide an even better 
possibility to study the spatial-temporal structure of the field QHorbury, 2000[ 
[Osman and Horbury, 2007}|Osman and Horbury, 2009a||Osman and Horbury, 2009bD .| 

The principal tool generally used to describe and measure field fluctuations is 
the two point correlation function. Our goal is to determine the complete three 
dimensional (3D) correlation functions of small scale IMF fluctuations. Previous 
studies only determined sections of the complete 3D function, used sometimes 
questionable symmetry assumptions and/or used relatively small data samples. 
In this paper we analyze more than 200 hours of data, deal with three axis 
anisotropy and reconstruct the complete 3D correlation function from its mea- 
sured sections. Using the correlation function we arc able to calculate for example 
three axes correlation ellipsoids, anisotropy and correlation length ratios. Our 
results can be used to verify theoretical predictions, to investigate turbulence 
and intermittency as well as the transport properties of the solar wind. 

In Section [5] we present the observations, in Section [3] the methods are de- 
scribed and finally we summarize our results in Sectional 
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Figure 1. The cylindrical projection of Cluster spacecraft orbits during our observations in 
GSE system, where X and R = \J^qse ^GSE ^'-^ measured in REarth units. A model 
magnetopause (left curve) and bow shock (right curve) arc also shown for reference. 

Table 1. The intervals studied. During these intervals the Cluster 
spacecraft dwelt inside the solar wind. 
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2. Observations 

For this study we have used Is averaged time resolution Cluster II Flux Gate 
Magnetometer (FGM) data ( [Balogh et ai, 2001[ ). The spacecraft separation was 
in the range of several hundreds to several thousands of kms during the obser- 
vations. The elongated ellipse orbit of Cluster spacecraft samples the solar wind 
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Figure 2. An example of the measured sections of the correlation function. Ci is the i*'' 
Cluster spacecraft. The red line (C1:C4) for example shows the correlation of the magnetic 
data measured on Clustcrl and Clustor4. 

in February, March and April and allows interplanetary magnetic field (IMF) 
measurements (Fig. [T]) during those months. One interval was chosen arbitrarily 
for each year between 2003 and 2008, in which the measurements had no data gap 
(See Table[T|) . The durations of the intervals are between 50 min and 28 hours. We 
also used Cluster Ion Spectrometer Hot Ion Analyzer (CIS/HIA) spin averaged 
resolution data ( |Reme et ai, 200"T| ) on Cluster SCI to determine the solar wind 
speed. We used this parameter to transform the data into the solar wind reference 
frame, which was the first step of the correlation calculation. 



3. Multi-spacecraft methods 

The purpose of this study is to determine full 3D correlation functions. Using 
the Is averaged FGM measurement we calculated one dimensional time-lagged 
auto- and cross-correlation functions of the magnetic fluctuations (Cluster Cl- 
Cl, C1-C2, C1-C3, C1-C4, See Fig. [5]) in every ten minutes of the measurement 
intervals of Table [TJ These are one dimensional sections of the 3D correlation 
function. After the calculation of the sections we fitted a single 3D test function 
to all of the points of all four sections. 

The correlation lengths are quite long in the solar wind but the separation 
of the Cluster fleet was enough to measure significant variations, therefore to 
calculate the correlations. Every section was calculated in a 120s long interval, 
i.e. the time lags was less then 60 seconds. Adaptive correlation calculation was 
used to decrease the calculation time: 

1. First a draft of the correlation function was calculated by 5 s long leapfrog 
steps. 

2. The maximum of the correlation functions was determinated. 

3. Finally the correlation function was calculated 10 s before and after the 
maximum using Is resolution. 

This method spared computational resources and focused on the most important 
(most rapidly changing) part of the correlation function. 
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Figure 3. Illustration of our method: The four Cluster spacecraft sampling a volume of the 
solar wind. We compute the correlations along ri-r,i, and fit a 3D correlation function to the 
resulting data. 



Using the Taylor's or "frozen in field" hypothesis ( [Taylor, 1938[ ), these one 
dimensional, time dependent measurements can be translated to samples of the 
3D spatial correlation function along a multitude of r„ space vectors (Fig. [3|) in 
solar wind frame. 

The Taylor's hypothesis states that on a time scale short compared to the 
characteristic timescale on which the fluctuations vary in the solar wind frame, 
the variations measured by the spacecraft are due to the convection, i.e. approx- 
imately static fluctuations are convected past the spacecraft. This is true only if 
the convection speed is much greater than any other speed influencing the mea- 
surements. Since both the velocity of the spacecraft and the wave propagation 
speed (roughly the Alfven speed Va) are quite low compared to the solar wind 
speed (Vsw) at 1 AU, the time lag (t) recorded by a spacecraft can be converted 
to a spatial position as x = V^^t ( [Taylor, 1938[ ). 

The Taylor's hypothesis however has a modified form for multi-spacecraft 
case. In this case the distance of two measurement positions in the solar wind 
frame can be quite low (even zero), thus the finite wave propagation speed 
becomes significant ( [Horbury, 2000[ ). In this case the condition for the Taylor's 
hypothesis can be written as |ri_2 — 'Vswt\ ^ IVa^I where ri^2 is the distance 
of the two spacecraft. The condition is fulfilled for the vast majority of our 
measurement points. For the few cases, when it is not, we have only a single 
corrupted point out of the 160 we use to fit a correlation function. (|ri.i — Vsmtl 
could only be small for one point for fixed ri.; vectors.) This single point has no 
significant influence on the fit. 

We can now fit a 3D test function (an Ansatz) C(x,y,z) to these samples 
(Fig. Uj) . The quality of the fit is measured by the sum of the squared deviations 
of our function and the measured values. We found that the function 



C {x,y,z) = exp 



f2 /2 /2 

X y z 
6^ q2 
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Figure 4. A two dimensional section of the 3D correlation function that we used to fit the 
data. The correlation values are given on the vertical axis and also by a color code with respect 
to the spacecraft separation in the X and Y directions 

fits very well in most of the time intervals we have studied. (Here a, b and c are 
the semi-principal axes of the correlation ellipsoids; x', y' and z' are transformed 
into the principal coordinate system.) 

Random processes can usually be described by an exponential correlation 

kl 

function C(r) oc e~~ , where A is the correlation length. Actually the correla- 
tion length is defined by this equation. (When the correlation function does not 
show exponential decay, the correlation length is infinite; these are the so-called 
strongly correlated systems, occurring usually only in very special situations 
such as phase transitions, critical points and alike.) 

When we set the 7 parameter to be 7 = 1/2 in our correlation function 
C {x,y, z), then that is only a 3D version of the usual equation, with three 
independent correlation lengths. Unfortunately if we fix 7 to a certain value, 
the quality of the fit is significantly worsened, so we were forced to use 7 as an 
independent parameter. 

When 7^1/2 the correlation function is a stretched-exponential. (This kind 
of correlation is often found in systems with broad distribution of correlation 
scales.) In this case the semi principal axes of the correlation ellipsoids are not 
exactly the classical correlation lengths, but their ratios serve as a very good 
measure of the anisotropy. 

We have seven independent parameters in a transcendent function (a, b, c, 7, 
and the three independent parameters of the rotation matrix) so we encountered 
serious problems during our analysis: 

• The Ansatz is a transcendent function so linear regression cannot be used. 

• The exponent is unknown so it cannot be solved for the logarithm. 

• The rotation matrix has 9 parameters but only 3 of them are independent. 

So our task was to fit a transcendent function with 9+3+1 parameters and 6 
constraints. 

We solved this task by using the least squares method in a special way. We 
re-defined the Ansatz by using a special parameter set taking into account the 
constraints. The new parameter set eliminates the usage of the transformation 
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Figure 5. Examples of the correlation ellipsoids described in Sec. [S] The ellipsoids are the 
" equipotential" surfaces for a fixed {C (x,y,z) = 0.99) value. The units are arbitrary but the 
size of the object is proportional. The correlation function is clearly essentially anisotropic. 



matrices in the parameterization of the correlation elhpsoids: 



Then we let the system move in the 7 dimensional parameter space on a surface 
defined by the least squares toward the minimum of the surface. Some random 
step was forced to avoid getting the local minima instead of the real minimum. 
Adaptive step size was used to spare resources and focus on the important region. 
Applying this method with these parameters the code finds stable minimum most 
of the time. The minimum does not depend on the choice of initial parameters. 
The final value of the error sum is 3-5 orders of magnitude smaller than its initial 
value. This value shows how well C{x,y,z) fits to the 160-160 measurement 
points of each evaluations. The sum is usually in the 10^^ — 10^^ range, which 
means that the error for a single measurement point is approximately 1— 3x 10~^. 
Since it is really difficult here to estimate the error from basic principles (among 
other things it depends on the relative positions of the spacecraft and the physical 
properties of the plasma) we used this error sum of the least squares fit to 
estimate the overall error of our method. 

Though the new parameters are extremely useful in solving the fitting prob- 
lem, it is not easy to find the physical meaning of them. But after the fit we 
can reconstruct the physical parameters for further analysis. From our special 
parameters we compute the magnitudes and directions of the semi principal axes 
of the correlation ellipsoids for C {x, y, z) = 0.99. 

4. Discussion and conclusion 

A lot of ellipsoids were visualized given by the algorithm described in Sec. |31 
These ellipsoids show significant variations in size as well as shape. Most of them 
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Figure 6. The distribution of the 7 exponent of the stretched-exponential correlation function 
(see Sec.[3]l, with its mean value and standard deviation. N is the number of evaluations with 
7 values falling in the corresponding 0.025 wide bin. 




Figure 7. Ratios (c:b, c:a, b:a) of the semi-principal axes of the correlation ellipsoids (i.e. 
ratios of the anisotropic correlation lengths). Samples from year 2003 to 2008. Black line: 
ratio=l. Purple line: fit for the first population. Orange line: fit for the second population. 



show strong three axes anisotropy. This means that the correlation functions and 
so the fluctuations of the solar wind plasma are basically anisotropic (Fig.[S]) they 
have neither spherical nor axial symmetry. 

The 7 exponent of the stretched-exponential correlation function shows a 
broad distribution. Its mean value and standard deviation are 0.33 and 0.12, 
respectively (Fig. IH]). 

We found that small scale IMF fluctuations are highly inhomogeneous: the 
correlation length varies over almost six orders of magnitude. The IMF turbu- 
lence shows significant anisotropy with two distinct populations. Most of the 
time the ratio of the 3 principal axes of the correlation ellipsoid, i.e. the ratio 
of the correlation lengths, is roughly 1:2.2:6, in the remaining time the major 
axes can be much larger (Fig. [T]). In this second population the ratio of axes 
shows high variations but the ratio 1:2.2:60000 has slightly higher probability. 
This means extremely strong correlation in a given direction. We found favoured 
directions in the orientation of the correlation ellipse (Fig. [5]). The major axes 
favor the Sun-Earth direction. The minor axes favour the direction perpendicular 
to the ecliptic plane. In most of the time, the correlation lengths computed using 
Bz data (GSE coordinates) are significantly smaller than those computed using 
Bx or By. The cause of this could be the relative lack of slab turbulence. 
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Figure 8. The directions of the three principal axes of the correlation ellipsoids. GSE coordi- 
nates, (a) X-Y, (b) X-Y and (c) Y-Z projections. Computed from data measured on 04/03/2005. 
Red: minor axes, blue: major axes. 



Acknowledgement The work of G. Facsko was supported by the OTKA grant K75640 of 
the Hungarian Scientific Research Fund. The authors thank Iannis Dandouras (CIS PI) and 
the Cluster Active Archive for providing CIS HIA solar wind measurements and Anna-Maria 
Vi'gh for improving the English of the paper. 



References 

Balogh, A., Carr, CM., Acuiia, M.H., Dunlop, M.W., Seek, T.J., Brown, P., Fornagon, K., 
Georgescu, E., Glassmeier, K., Harris, J., Musmann, G., Oddy, T., Schwingenschuh, K.: 
2001, The Cluster Magnetic Field Investigation: overview of in-flight performance and 
initial results. Annales Geophysicae 19, 1207-1217. 

Bieber, J.W., Wanner, W., Matthaeus, W.H.: 1996, Dominant two-dimensional solar wind 
turbulence with implications for cosmic ray transport. J. Geophys. Res. 101, 2511 — 2522. 
doi:10.1029/95JA02588 

Bigazzi, A., Biferale, L., Gama, S.M.A., Velli, M.: 2006, Small-Scale Anisotropy and In- 
termittence in High- and Low-Latitude Solar Wind. Astrophys. .]. 638, 499-507. 
doi: 10.1086/498665 

Dasso, S., Milano, L.J., Matthaeus, W.H., Smith, C.W.: 2005, Anisotropy in Fast and Slow 
Solar Wind Fluctuations. Astrophys. J. Lett. 635, L181-L184. doi: 10.1086/499559 

Horbury, T.S.: 2000, Cluster II Analysis of Turbulence using Correlation Functions. In: 
R. A. Harris (ed.) Gluster-II Workshop Multiscale / Multipoint Plasma Measurements, 
ESA Special Publication 449, 89. 

Matthaeus, W.H., Goldstein, M.L., Roberts, D.A.: 1990, Evidence for the presence of quasi- 
two-dimensional nearly incompressible fluctuations in the solar wind. J. Geophys. Res. 
95, 20673-20683. doi 10.1029/JA095iA12p20673 

Matthaeus, W.H., Dasso, S., Weygand, J.M., Milano, L.J., Smith, C.W., Kivelson, M.G.: 2005, 
Spatial Correlation of Solar-Wind Turbulence from Two-Point Measurements. Physical 
Review Letters 95(23), 231101. doi: 10.1103/PhysRevLett.95. 231101 

Nemeth, Z., Erdos, G., Balogh, A.: 2000, Fluctuations of the Hcliospheric Magnetic 
Field and the Anisotropy of Pickup Ions. Geophys. Res. Lett. 27, 2793-2796. 
doi:10.1029/2000GL000068 

Nemeth, Z., Facsko, G., Erdos, G., Balogh, A.: 2006, Plasma mixing as a cause of 
solar wind magnetic field variations. Advances in Space Research 37, 467-472. 
doi:10.1016/J.asr.20 05.11.012| 

Osman, K.T., Horbury^ T.S.: 2007, Multispacecraft Measurement of Anisotropic Corre- 
lation Functions in Solar Wind Turbulence. Astrophys. J. Lett. 654, L103-L106. 
doi: 10.1086/510906 

Osman, K.T., Horbury, T.S.: 2009a, Multi-spacecraft measurement of anisotropic power levels 
and scaling in solar wind turbulence. Annales Geophysicae 27, 3019-3025. 

Osman, K.T., Horbury, T.S.: 2009b, Quantitative estimates of the slab and 2-D power in 
solar wind turbulence u sing multispacecraft dat a. Journal of Geophysical Research (Space 
Physics) 114, 6103. doi |10.1029/2008JA014036| 



SOLA: sola-s-10-00092-4-tx.tex; 17 June 2010; 1:54; p. 9 



Nemeth, Facsko, and Lucek 



Remo, H., Aoustin, C, Bosqucd, J.M., Dandouras, I., Lavraud, B., Sauvaud, J. A., Barthe, A., 
Bouyssou, J., Camus, T., Coeur-Joly, O., Cros, A., Cuvilo, J., Ducay, F., Garbarowitz, 
Y., Medale, J.L., Penou, E., Perrier, H., Romefort, D., Rouzaud, J., Vallat, C, Alcayde, 

D. , Jacquey, C, Mazelle, C, D'Uston, C, Mobius, E., Kistler, L.M., Crocker, K., Granoff, 
M., Mouikis, C, Popecki, M., Vosbury, M., Klecker, B., Hovestadt, D., Kucharek, H., 
Kuenneth, E., Paschmann, G., Scholer, M., Sckopke, N., Seidenschwang, E., Carlson, C.W., 
Curtis, D.W., Ingraham, C, Lin, R.P., McFaddcn, J. P., Parks, G.K., Phan, T., Formisano, 
v., Amata, E., Bavassano-Cattaneo, M.B., Baldctti, P., Bruno, R., Chionchio, G., di Lellis, 

A. , Marcucci, M.F., Pallocchia, G., Korth, A., Daly, P.W., Graeve, B., Rosenbaucr, H., 
Vasyliunas, V., McCarthy, M., Wilber, M., Eliasson, L., Lundin, R., Olsen, S., Shelley, 

E. G., Fuselier, S., Ghielmetti, A.G., Lennartsson, W., Escoubet, CP., Balsiger, H., Friedel, 
R., Cao, J.B., Kovrazhkin, R.A., Papamastorakis, I., Pellat, R., Scudder, J., Sonnerup, 

B. : 2001, First multispacecraft ion measurements in and near the Earth's magnetosphere 
with the identical Cluster ion spectrometry (CIS) experiment. Annales Geophysicae 19, 
1303-1354. 

Richardson, J.D., Paularcna, K.I.: 2001, Plas ma and magnetic field correlations in the solar 
wind. J. Geophys. Res. 106, 239-252. doi:'lD.1029/2000JA000071 

Shcbalin, J.V., Matthaeus, W.H., Montgomery, D.: 1983, Anisotropy in MHD turbulence due 
to a mean magnetic field. Journal of Plasma Physics 29, 525-547. 

Taylor, G.I.: 1938, The Spectrum of Turbulence. Royal Society of London Proceedings Series 
A 164, 476-490. 



SOLA: sola-s-10-00092-4-tx.tex; 17 June 2010; 1:54; p. 10 



